The ^S+^S asymptote of Sr2 studied by Fourier-transform spectroscopy 
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An experimental study of the long range behavior of the ground state X^E^ of Sr2 is performed 
by high resolution spectroscopy of asymptotic vibrational levels and the use of available photoassoci- 
ation data. Ground state levels as high as v" = 60 (outer turning point at 23 A and 0.1 cm"'^ below 
the asymptote) could be observed by Fourier-transform spectroscopy of fluorescence progressions 
induced by single frequency laser excitation of the v' = 4, J' = 9 rovibrational level of the state 
2^EJ. a precise value of the scattering length for the isotopologue ®*Sr2 is derived and transferred 
to all other isotopic combinations by mass scaling with the given potential. The derived potential 
together with already published information about the state 2^Ej directs to promising optical paths 
for producing cold molecules in the electronic ground state from an ultracold ensemble of Sr atoms. 

PACS numbers: 34.20.Cf, 31.50.Bc, 33.20.Kf, 33.20.Vq 



I. INTRODUCTION 

Our spectroscopic work on this molecule is motivated 
by the current high interest in ultracold ensembles of 
Strontium atoms [1-4], vi^hich could be a candidate for 
an optical frequency standard |5,;J1[, and also by inter- 
est in ultracold Sr2 molecules [3, 1^. Though for all 
these experiments reliable knowledge of collision prop- 
erties for Sr atoms like scattering lengths would be of 
advantage, there was no sufficiently precise ground state 
potential available from ab initio calculations or exper- 
imental work, from which these quantities could be de- 
rived by directly solving the radial Schrodinger equation. 

In [lOf we reported on the optical spectrum of the sys- 
tem 2^S+ — ^^^t ^^^^ more than 10300 transitions 
which involve the term energies of nearly 60% of the ex- 
isting rovibrational levels of the ground state for the main 
isotopologue ^*Sr2 and additionally many levels for the 
isotopologues ^^Sr^^Sr and ^''Sr^^Sr. We reached vibra- 
tional levels up to v" — 49 and rotational levels up to 
J" = 224. Using this data set we had to correct the ro- 
tational assignment given by [ll[ by four units and were 
able to construct a potential energy curve (PEC) which 
is reliable in the range from 4 to 11 A, covering an energy 
region from the bottom up to 9 cm~^ below the dissocia- 
tion asymptote. With the help of theoretically calculated 
long range coefficients Gq and Cg (l^ . [T3| and at that time 
available node positions of the scattering wave function at 
large internuclear separation for a kinetic energy of a few 
microkelvins [l^ and 2 mK [ij from photoassociation, we 
estimated a complete set of scattering lengths for all com- 
binations of naturally abundant Sr isotopes. Using our 
derived potentials to calculate Franck-Condon factors we 
proposed spectroscopic schemes to measure the highest 
existing vibrational levels of the ground state. Meanwhile 
new work on photoassociation 15] was performed and the 
binding energy of the last vibrational level v" = 62 for 
**Sr2 was determined and for ^^Sr two groups reached 
Bose-Einstein condensation [l^, [l3| . 

In this work we report on successful experiments in the 
proposed direction and derive ground state levels close to 



the asymptote with which we extend the precisely known 
potential region significantly into the long range regime. 
With these data the experimental knowledge about the 
long range coefficients Ce, Cg and Cio and the precision 
in the scattering lengths will be improved. 

The paper is organized as follows: Section HIl describes 
the experimental methods. Sec. IIIII shows the obtained 
data set and discusses the measurement uncertainties. 
Sec. IIVI presents the resulting ground state potential and 
scattering lengths. Sec. |V] gives conclusions and a short 
outlook. 



II. EXPERIMENT 

The apparatus is the same as in (loj , the experimental 
procedure is an extension of the earlier one. Strontium 
is filled into a stainless steel heatpipe, which is heated to 
a temperature of about 1220 K under 20 mbar of argon 
as a buffer gas. To keep the optical path free from con- 
densing strontium crystals growing at both ends of the 
heatpipe, the oven has to be moved every two hours to 
bring this Sr back to the heated zone. This means, that 
the heated cell is not working perfectly as a heatpipe 
where the condensed substance flows back to the heated 
region. 

The molecules are excited to the v' = 4 level of the 
state 2^1]+ by light from a Coherent CR 699 ring dye 
laser operated with Rhodamine 6G at frequencies close 
to 17635 cm~^ (depending on the rotational quantum 
number J). The power of the dye laser was about 70 mW 
in front of the heatpipe. The fluorescence light emitted 
antiparallel to the beam direction of the laser is imaged 
into a Bruker IFS 120 HR Fourier-transform spectrome- 
ter. 

To measure the highest vibrational levels of the ground 
state, the vibrational band of the excited state with 
v' = 4 and the rotational levels with J' — 25, J' — 21, 
J' = 17, J' = 13, and J' = 9 were selected. Low J are very 
important for observing asymptotic levels because high J 
are not existing for these v. For measuring the progres- 
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sion starting from J' = 25 a resolution of 0.05 cm"-'^ was 
used, but because of the small spacings between the lines 
with v" > 50 it was expected that for the smaller J' the 
P, R doublets of the highest vibrational levels would not 
be completely resolved with this resolution. Therefore, 
the resolution was increased for the remaining measure- 
ments to 0.02 cm~^. For increasing the sensitivity of 
the measurements without saturating the amplifiers in 
the detection unit of the spectrometer selected optical 
filters were used. For the progressions of the four larger 
J' (25 to 13) a holographic notch filter and a OG 590 
color glass filter were inserted to suppress the internal 
HeNe laser of the spectrometer and the stray light of the 
dye laser, respectively. For the three progressions (J' — 
25, 21, and 17) up to 600 scans of the spectrometer were 
averaged and for J' = 13 about 800 scans were necessary 
to identify safely lines to the highest existing levels with 
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= 59 and 60. To be able to observe 



= 60 also 



for J' — 9, where the intensity is already significantly 
lower because of the smaller degeneracy of that level, the 
experiment had to be further optimized. With a single 
bandpass fiUer with 10 nm FWHM (full width at half 
maximum) around 600 nm, where the expected fluores- 
cence for v" > 30 should appear, and a laser power of 
more than 200 mW in front of the heatpipe it was pos- 
sible to identify two more rotational levels for v" = 60 
(J" — 8 and J" — 10) by averaging 600 scans. 

The method, mentioned in [lOj and successfully used 
for Ca.2 in jl^] , of excitation spectroscopy by tuning the 
laser over the transition frequencies to the desired levels 
and filtering the detected fluorescence with a monochro- 
mator, was not used for observing even higher v", i.e. 
61 and 62, because for these levels, which exist only for 
J" < 8, the line separations of the P, R doublets become 
already smaller than the Doppler widths of the individ- 
ual lines. Beam experiments would be needed to obtain 
sufficient resolution for this purpose. 



III. DATA SET 

The data set published in [l3] consisted of 4673 dif- 
ferent rovibrational levels for the ground state, contain- 
ing 3163 of the ca. 5400 existing ground state levels of 
the main isotopologue ®®Sr2. The current ground state 
data set is extended to 6077 different rovibrational lev- 
els, of which 4159 belong to the isotopologue *^Sr2 and 
1301 and 617 belong to the isotopologues ®^Sr^®Sr and 
^^Sr^^Sr, respectively. The complete data field is shown 
in figure [T] 

The assignment of the 15700 lines (5400 new since pub- 
lication of 10]) of the current data set was done using a 
partially automated software as already described in [13] 
(in a prior version in [l^ ) . This software uses the differ- 
ences between individual lines of a selected progression 
to automatically assign them to a molecule and an elec- 
tronic transition and to quantum numbers v', J', v" and 
J", provided that the data of the ground and as much as 




FIG. 1: (Color online) Overview of the observed energy levels 
of the ground state X^Ei", including all measurements of ref. 



available of the excited states of the molecules are input 
to the program. Additional lines which cannot be auto- 
matically identified to belong to the progression, because 
they extend significantly beyond the precisely known re- 
gion of the ground state potential, are individually se- 
lected by the user before the program automatically as- 
signs the best fitting quantum numbers. 

In [13] the uncertainties for the differences of the line 
frequencies were estimated to be 0.005 cm"-'^ in most 
cases. This is one tenth of the applied resolution of 0.05 
cm~^. We keep this uncertainty for the lines which are 
taken from the earlier data set and also the lines newly 
measured using the same resolution. For the new lines 
measured with the higher resolution of 0.02 cm"'^ we esti- 
mate a smaller uncertainty of 0.002 cm~^ for the stronger 
lines (S/N ratio larger than 5). For these lines a signifi- 
cant contribution of any Doppler shift can be excluded, 
since the excitation frequency was optimized on the cen- 
ter of the Doppler profile and thus the velocity compo- 
nents of the excited molecules parallel to the direction 
of observation will be very small. The specified uncer- 
tainty of the Fourier spectrometer for frequency differ- 
ences within a spectrum is 0.001 cm^^, still below the 
assumed experimental uncertainty. For weak lines and 
lines belonging to accidentally overlapping excitations an 
uncertainty of 0.004 cm~^ is assumed. 

The complete list of assigned lines can be found in the 
additional online material. 



IV. RESULTS 

For the ground state a potential energy curve is cal- 
culated using the same type of representation as in [loj 
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(details can be found in [20] )■ The central part of the 
potential {Ri < R < Ra) is described by 



(1) 



J>1 



with the nonlinear mapping function 



R — R„ 



R + bR„i 

The inner repulsive wall (i? < -Ri) is given by 

^^W = ^ + |r' 



(2) 



(3) 



while the long range part (i? > Ra) is represented by 



VaiR) - t/oo - 



i?6 



i?8 



(4) 



The coefficients Uj arc the main fitting parameters, 
while the parameter T„i is adjusted to obtain at the point 
Ra a continuous connection between the central and the 
long range part of the potential. Rm in eqn. ([2]) is chosen 
close to the minimum of the potential. The parameters 
b and Ri are manually chosen to optimize the fitting re- 
sult with a small number of parameters Uj. The real 
parameter n [2l| is an additional fitting parameter, while 
the parameters A and B are adjusted for a continuously 
differentiable connection at the point Ri. 

The radius Ra should be close to the LeRoy radius 

Rb{AB) = 



22l |. where ta and tb are 



the radial coordinates of the outermost electrons of atom 
A and atom S, here = = rsr- The resulting LeRoy 
radius Rb is 9.7 A or 10.5 A for the X state of Sr2, if (r|^) 



is taken from 23] or from [24|, respectively. The larger 
value Ra — 10.5 A was chosen for the final potential 
evaluation. This choice should justify the neglect of a 
formal term for an exchange energy in the long range 
part according to cq. 

The residuals of the potential fit can be found as his- 
tograms in figure [2] compared to Gaussian profile fits. 
The plot for the main isotopologue **Sr2 shows a sym- 
metric profile with a FWHM of 0.8, while the shapes for 
the other two isotopologues are significantly wider and 
are not as symmetric as for the main isotopologue. This 
can be attributed to the much lower number of data for 
the isotopologues ^^Sr^^Sr and ^^Sr^^Sr and on average 
significantly weaker line intensities, which lead to a less 
precise determination of the frequencies. 

The derived potential coefficients are given in table ID 
The normalized standard deviation of this potential fit 
is (7 = 0.80 which indicates that the applied measure- 
ment uncertainties are as a whole not underestimated. It 
correlates to the width of the histogram in fig. [5] which 
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FIG. 2: (Color online) The residuals of the X^Eg potential 
fit divided by the specific uncertainties represented by his- 
tograms compared to Gaussian fits. 



also gives a normalized FWHM of 0.8 for the main iso- 
topologue. The values of the dissociation energy De and 
the equilibrium internuclear distance Re depend on the 
chosen potential description, namely on the number of 
potential coefficients Oj, on the parameter b and on the 
connection radii Ri and Ra- The given uncertainties are 
Icr standard deviations derived by comparing 35 differ- 
ent potential versions with different parameter settings 
which all had almost the same a values. Compared to 
De the value for Do (dissociation energy with respect to 
(v=0,J=0) for the main isotopologue *^Sr2) is less model 
dependent. 

For the first potential fits the long range coefficients 
Cg, Cg, and Ciq were treated as free parameters. De- 
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TABLE I: The potential coefficients and derived molecular 
constants for the ground state X^E^. The potential energy is 
calculated with respect to the asymptote, the dispersion co- 
efficients are the recommended values from the third column 
of table El 



ai 


-6.50x10"^ 


cm~i 


a2 


1.5939056x10'' 


cm"'^ 


as 


-2.9646778x10'' 


cm"'^ 


34 


-6.269777x10^ 


cm"'^ 


as 


4.4952358x10" 


cm"'^ 


ae 


8.709016x10^ 


cm"'^ 


a? 


-1.0054929x10^ 


cm"'^ 


as 


5.94784152x10^ 


cm"'^ 


ag 


-9.95239126x10^ 


cm"'^ 


aio 


-1.14496717x10^^ 


cm"'^ 


an 


4.606463055x10^ 


cm"'^ 


ai2 


3.74666573x10^ 


cm"'^ 


ai3 


-5.439157146x10* 


cm"'^ 


ai4 


9.364833940x10** 


cm"'^ 


ais 


1.387879748x10^ 


cm"'^ 


ai6 


-8.4009054730x10^ 


cm"'^ 


ai7 


1.5781752106x10^° 


cm"'^ 


ai8 


-1.5721037673x101" 


cm"'^ 


ai9 


8.376043061x10^ 


cm"'^ 


a.20 


-1.88984880x10^ 


cm"'^ 


b 




-0.17 


Rm 


4.6719018 A 




-1081.6384 


cm 1 


Ri 


3.963 A 


n 


12.362 


A 


-1.3328825x10^ 


cm ^ 


B 


3.321662099x10^° cm'^A" 


Ra 


10.5 A 


Ce 


1.525x10^ cm^^A^ 


Cg 


5.159x10** cm-^A* 


Cio 


1.91x10^° cm 


-1^10 


derived constants: 






De 


1081.64(2) 


cm^'^ 



Do 1061.58(1) cm-i 
Rs 4.6720(1) A 



spite the precision of the energies of the levels v" — 62, 
J" — and J" — 2 from [lE\ incorporated in the fit, the 
set of coefficients is still significantly correlated. While 
the Cg coefficient is quite well fixed the coefficients Cs or 
Cio can be varied over a larger range of values without 
a significant reduction in the quality of the description 
of the data set, if all other coefficients are adjusted ac- 
cordingly. To estimate uncertainties for the long range 
coefficients different potential fits were done where each 
time the Cg coefficient was held fixed at a different value 
within its range as given in the second column of table 



im The results shown for the other Ci coefficients are av- 
eraged values from the series of fits and their uncertainty 
bounds represent an increase of the sum of the weighted 
squared residuals by 2%, which is a quite safe limit, since 
only 113 of the 6077 energy levels included in the fit are 
directly affected by this change. 

Because the derived uncertainty for the Cs coefficient 
is significantly larger than the uncertainty claimed for 
the theoretically calculated coefficient by [l^ (see the 
fifth column of table [11]) and also larger than the dif- 
ference between the two available theoretical Cg coeffi- 
cients [13, [13, additional potential fits were done, where 
each time the Cs coefficient was held fixed at a value 
taken from [l^l or [l^ or the average of both. Such a 
procedure was also tried with the theoretical Ce coeffi- 
cients, but those resulting potentials gave a significantly 
less convincing description of the data set and also the re- 
sulting coefficients Cg and Cio showed strong deviations 
to the theoretical values. Additional fits were done for 
different numbers of potential coefficients aj and differ- 
ent connection radii Ri and Ra ■ For a few fits Ra was set 
to values as large as 13 A, which led to significantly less 
well determined long range coefficients with higher cor- 
relation without obtaining a significant improvement in 
the description of the data. The resulting Ci coefficients 
of all the fits using either the Cg coefficient from or 
from ^13i] or its mean value were averaged to obtain the 
recommended values, which are listed in the third col- 
umn of the table [TTl The uncertainties are la standard 
deviations of the averaged values. 

Table |TT] compares the dispersion parameters derived 
with freely varied C; coefficients and the recommended 
values with the published theoretical long range coeffi- 
cients from recent years. As one can see, the values in 
the second column are in good agreement with all the 
other values, while the recommended values do deviate 
from the most precise theoretical coefficient set from [l^ 
more than one would expect considering the given un- 
certainties, there. This could be due to different reasons. 
The paper [3] used the lifetime of the atomic — 
transition from [l3|j which was calculated from the C3 
coefficient of the state 2^S]+ fitted to photoassociation 
data. In 14] line broadenings of the PA resonances are 
mentioned, which could indicate predissociation due to 
the coupling to the continuum of electronic states from 
lower asymptotes and could lead to a systematic shift of 
these asymptotic levels of the state 2^S+. Another rea- 
son could be that the presently derived uncertainties are 
underestimated because of a too low number of differ- 
ent potential descriptions used for averaging. The com- 
puting effort of testing the long range correlation by a 
Monte Carlo study would be enormous for the presently 
employed analytical type of potential description and the 
available hard and software. 

A pointwise representation of the potential derived 
from the parameters in table |T] and the analytical po- 
tential description for the best fitting values obtained by 
freely varying Ci coefficients can be found in the addi- 
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TABLE II: Comparison of the long range coefficients derived from fits of freely varying parameters with the recommended 
values applying a fixed Cg value and with theoretical long range coefficients published in recent years. For the details on 
evaluation procedure and uncertainty estimates see text. 





freely var. 


recommended 


[12] [25] [26] [27] 


Cs [lO^cm-iA^] 


1.527(33) 


1.525(5) 


1.566 1.495(3) 1.5481 1.528(94) 1.509(20) 


Cs [10«cm-iA«] 


5.1(15) 


(5.159) [28] 


5.201 5.118(4) 


Cio [10i°cm-iA"'] 


2.0(13) 


1.91(27) 


1.61 1.593 



tional online material. 

For the convenience of the reader we also made fits us- 
ing the MLR potential description from [2^ , which might 
allow for defining the potential by a lower number of pa- 
rameters. It has the mathematical form 



Vmlr{R) = D. 
with the long-range function 



ULRiRe) 



(5) 



Ulr{R) 



Cp, Cs 



C 



10 



i?6 ' i?8 ' i^lO 

and the exponential function 



(6) 



N 



13 = y;-f iR)f3^ y;^f (i?)] ^ (r)^ , (7) 

where 



/3oo = In 



2D, 



ULR{Re)_ 

Three mapping functions y'^ are defined as 



(8) 



Vb 



R-Rl 
R" + Ri • 



(9) 



where the placeholder a can either be 'e' for the equi- 
librium radius Re or Ve/' for Rref and the exponent b 
can either he p or q. Rref is an expansion point at an 
internuclear distance larger than Rg. The f3i, Ci, Re, and 
the dissociation energy De are the fit parameters, while 
the parameter Rref and the exponents p and q are man- 
ually chosen to optimize the potential description with 
a minimum number of coefficients and to get the proper 
long range behavior. 

The potential coefficients for the MLR description are 
given in table IIIIl the normalized standard deviation is 
a — 0.80 as for the description given in table HI Thus 
there is no significant difference in the quality of the de- 
scription of the data set. The clear advantage of the 



TABLE III: The coefficients of the MLR potential description 
for the ground state X^E^. The potential energy is calculated 
with respect to the potential minimum, the dispersion coeffi- 
cients were taken from table HIl 



A) 


-1.458760592449 


I3i 


-0.031511799120 




-0.912033020316 




-0.033913017557 




-0.465124882150 


P5 


0.287008644909 


/36 


-0.585285048765 




2.490073196261 




3.065232204927 




-8.521057394373 


/3io 


-0.771989570637 


/3ii 


17.679758812808 


/3l2 


-11.143746759192 


I3l3 


-1.704727367184 


Ce 


1.525x10'^ cm-^A*^ 


Cs 


5.159x10* cm-^A* 


Cio 


1.910x10^" cm-^A^" 


De 


1081.6352 cm"^ 


Re 


4.6719032 A 


Rref 


5.5 A 


P 


5 


q 


3 



potential description from table IIIIl compared to the po- 
tential description from table H] is the lower number of 
potential coefficients necessary (14 /3i instead of 20 a^) 
and the fact that it is indefinitely continuously differen- 
tiable at all points which is not the case for the poten- 
tial of table U at the points Ri and Ra- Disadvantages 
are the much more complex formula with three mapping 
functions instead of one and a much stronger correlation 
of the long range coefficients Ci with the potential co- 
efficients f3i (see eqn ([5])). For the potential description 
from table U the correlation between the Ci and the is 
only due to the integration of the wave functions over the 
full potential region. Thus we applied the recommended 
Ci values from table |ll] for the fit of the MLR potential. 
The dissociation energy obtained in the MLR approach 
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agrees very well to the derived value in table ID 

Table |TV] shows the derived scattering lengths of the 
different isotopic compositions and compares them to 
other recent publications. The values in the second col- 
umn were calculated from the potentials which were ap- 
plied for the estimation of the uncertainty of the long 
range coefficients as given in the second column of ta- 
ble ini The values in the third column of table IIVI are 
averages from different potentials which were obtained 
using Cg coefficients from 12 1 and [l^ and resulted in 
the recommended values in the third column of table |lll 
ft turned out that the influences of different connection 
radii Ri or Ra or numbers of parameters Oj , which causes 
mainly a change in the central part of the potential, has a 
much larger influence on the scattering leng ths than the 
question if the Cs should be taken from [13 or from p^ . 
All of the potentials used reproduce the measured data 
nearly equally well and also the node positions of the 
scattering wave functions from photoassociation [ll, [3] ■ 
The uncertainties given in the second and third column 
of table IIVI are la standard deviations in units of the 
last digit shown. The uncertainties in the fourth column 
are taken from fis'l . The last column provides the results 
from our earlier analysis in [lo| for comparison. The very 
precise analysis of the present work is consistent with the 
values froni[10| and with those from the photoassocia- 
tion work [15j . The new determination yields about a 
factor of 2 smaller uncertainties than reported in . 



V. CONCLUSION 

After publishing a description for the ground state po- 
tential for the region from 4 to 11 A and the vibrational 
interval from v" = to 48 in [l^l, for the present pa- 
per we concentrated our work on the asymptotic levels 
of this state. By applying the same spectroscopic method 
but improving the detection scheme, increasing the signal 
averaging time and selecting appropriate excited rovibra- 
tional levels to reach the long range regime we were able 
to observe levels with v" up to 60, while v" — 62 would 
be the last existing one for low J". This allowed us to ex- 
tend the potential region up to an internuclear distance 
of 23 A from our data set. 

Including the term energies of the v" = 62, J" = 
and J" = 2 levels measured by photoassociation [Tsj, we 



were able to derive the long range coefficients Cq , Cs and 
CiQ. By taking the Cg coefficient from theory OIIIIj 
which has in our modeling a larger uncertainty than the 
theoretically calculated Cg coefficient, we derived an in- 
ternally consistent set of the remaining two coefficients 
Cg and Cio- This leads to a further reduction of the un- 
certainties in the calculation of the scattering lengths for 
all combinations of natural abundant isotopes of Stron- 
tium. We were also able to give a reliable value for the 
dissociation energy Dq of the X state. 

The combination of the data of both the state 2^E+ 
[To| . and the ground state X-^S+ can be useful for ex- 
periments producing ultracold molecules. The transition 
path successfully applied for this work to observe vibra- 
tional levels of the ground state close to the asymptote 
can be used in the reverse direction to transfer popula- 
tion of highly excited ground state levels to deeply bound 
vibrational levels of the X state within a single STIRAP 
step [13]. The asymptotic ground state levels, which can 
be populated e.g. by spontaneous decay after photoasso- 
ciation, show the best Franck-Condon overlap with the 
v' = 4, v' = 3 and v' = 9 levels of the state 2^1]+. While 
either the v' = 4 or the v' = 3 level give good overlap 
with almost every level of the X state, by excitation of 
the v' = 9 level almost 12% of the spontaneous decay 
reaches the v" = level. By the use of the 2^1]+ po- 
tential determined in [l3| together with the ground state 
potential from this work the transition frequencies from 
all ground state levels to the minimum region (v'=0 - 12) 
of the state 2^1]+ should be easily calculable with a preci- 
sion better than 0.01 cm~^. This accuracy should reduce 
substantially the time for finding STIRAP transitions. 

For the near future we plan to further investigate the 
excited states of the Sr2 molecule. We already started 
with the state 2^1]+ and we will look for the state l^Sj 
and a state of type for which we found one Q pro- 
gression during this work. 
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